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Abstract 

The critical values of the Yang- Yang functional corresponding to the vacuum states 
of the sine-Gordon QFT in the finite-volume are studied. Two major applications are 
discussed: (i) generalization of Fendley-Saleur-Zamolodchikov relations to arbitrary 
values of the sine-Gordon coupling constant, and (ii) connection problem for a certain 
two-parameter family of solutions of the Painleve III equation. 
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1 Introduction 

Throughout the past, a number of important facts about the quantum sine-Gordon model were 
discovered. Among them are elegant relations between the zero-point energy and the Painleve 
III transcendent. To describe them explicitly let us recall some elementary facts about a 
structure of the Hilbert space of the model, 

c= w ( \ (<9 ^ )2+A cos(0) ) ' (L1) 

in finite-size geometry with the spatial coordinate x compactified on a circle of a circumference 
R, with the periodic boundary conditions 

<f>(x + R,t) = 4>(x,t) . (1.2) 
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Due to the periodicity of the potential term A cos(0) in (11 .ip in 0, the space of states H splits 
into orthogonal subspaces Hk, characterized by the "quasi-momentum" k, 

0^0 + 2vr : | * fc ) -> e 2 ^ ifc | * fc ) (1.3) 

for | } e We call fc- vacuum the ground state of the finite-size system (11. ip in the sector 
1-ik and denote it by | v|/j s vac - ) y The corresponding energy will be denoted by E^. 

In general, the coupling constant in ( 11. ip should be restricted by the condition [3 2 g < 8n [1] 
and it is convenient to substitute (3 2 g for the "renormalized coupling" : 

B 2 

e = sg . (i.4) 

The value £ = 2 is special. For this coupling, the theory possesses M = 2 supersymmetry which 
is spontaneously broken, except the subspaces %k corresponding to k — ±~ [2]. In the sectors 
with unbroken supersymmetry the ground state energy is of course identically zero. In Ref. [3] 
Fendley and Saleur (see also related Ref. (3J) applied the general construction [5] to derive the 
remarkable relation 

R fdE k \ A dU(r) 

Here the variable r stands for the size of the system measured in the units of the correlation 
length (inverse soliton mass M), 

r = MR , (1.6) 
and U — U(t) is a particular solution to the Painleve III equation 

1 d/ t^) = 1 - sinh(2C/). (1.7) 



t dt V dt / 

This equation admits a one-parameter family of solutions regular at £ > (see e.g. [G]) called 
the Painleve III transcendents. The special solution in ( II. 5p is fixed by the following boundary 
conditions 

f-jlog(t) + 0(l) as t-0 
I o(l) as £ — > oo 

In the consequent work [7] Alyosha Zamolodchikov derived one more mysterious relation 

R />)Ek ^ - r l + l r dtt smh 2 U(t). (1.9) 



tt V d£ J Jf 1/4 8 2 

Below, we will refer to Eqs. (II. 5p . ( II. 9p as the FSZ relations. 

Relations similar to ( II. 5p and (II. 9p were also discovered in other models [1], [8], [9]. How- 
ever all the generalizations had limitations in the choice of coupling constants and sectors of 
the theories. The long-time consensus was the FSZ relations are due to the accidental sym- 
metry and do not possess any interesting generalizations for general values of £ and k. The 
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first serious sign that this may not be true came from the study of D = 4 supersymmetric 
gauge theories [TTHT5] . In these works a link was observed between certain Thermodynamic 
Bethe Ansatz (TBA) type integral equations and partial differential equations integrated by 
the inverse scattering methods. Some of the integral equations were in fact identical to the 
sine-Gordon TBA systems corresponding to £ 7^ 2 and k 7^ ±|. Inspired by this remarkable 
development, A. Zamolodchikov and the author found a classical integrable equation associated 
with the quantum sine-Gordon model for generic £ and k [TO]- It turned out to be the classical 
Modified Sinh-Gordon equation (MShG) 

d z d s 7] - e 2r > + p(z) p(z) e- 2 " = (1.10) 

with p(z) of the form 

p( z ) = z 2a - s 2a . 

Parameters a and s are real and positive, related to the sine-Gordon parameters £ (11. 4p and 
r = MR flU]) as follows 
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where, for future references, we use the notation 

r £ = ^%. (1.12) 

r(| + |) 

The MShG equation in general has no rotational symmetry. Instead, it has the discrete sym- 
metry z — > e^ z , z — > z. Solutions of the MShG equation (ll.lOp relevant to the problem 
respect this symmetry, are continuous at all finite nonzero z, and grow slower then the ex- 
ponential as \z\ 00. In other words, they are single-valued functions on a cone with the 
apex angle - including the zero of p(z) (see FigJT]). There is a one-parameter family of such 
solutions, characterized by the behavior at the apex: 77 — )- 2Z log \z\ +0(1) as \z\ — > 0, with real 
I E (— ~, I ) which turns out to be related to the quasi- momentum (jl.3p by 

I = 2\k\ - \ . (1.13) 

The MShG equation can be represented as a flatness condition for certain SL(2) connec- 
tion. In Ref. [10] it was shown that the monodromy from the apex to infinity corresponding 
to the above described solution is essentially the sine-Gordon Q-function, whose asymptotic 
expansions generate the vacuum eigenvalues of integrals of motions of the quantum theory. 

The original motivation for the present work was to incorporate the FSZ relations to the 
construction of Ref. [TU]. This problem is solved in Section[5]of this work. It turned out that the 
main player in the derivation of the generalized FSZ relations is a properly defined "on-shell" 
action for the MShG equation. Remarkably it can also be interpreted as a critical value of the 
Yang- Yang (YY) functional in the quantum sine-Gordon model. 

In the seminal work [16] the variational principle was applied to prove an existence of a 
solution to vacuum Bethe Ansatz (BA) equations for the XXZ spin chain. From that time the 



3 




Figure 1: The world sheet for the MShG equation f ll.lOj) . The dots A and B indicate positions 
of the apex and zero of p(z), respectively. At the minisuperspace limit (a — > oo, r is kept fixed) 
the world sheet shrinks to a single ray. 



functional whose extremum condition reproduce BA equations bears the Yang- Yang name. The 
YY-functional proves to be useful for computing norms of the Bethe states (see e.g. [Uj and 
references therein). Recently it attracts a great deal of attention in the context of the relation 
between supersymmetric gauge theories and quantum mechanical integrable systems [T8l - f20] . 
However, the role of the YY-functional in 2D QFT seems to be undervalued. To the best of 
my knowledge, it was never defined in intrinsic terms of integrable QFT. Nevertheless, there 
is a brute-force approach for the calculation of critical values of the YY-functional in the sine- 
Gordon QFT. It is based on the discretization of the theory, i.e., reducing it to the system 
with finite number of degrees of freedom, which then can be solved by standard methods of 
BA. Although this formal approach does not clarify the meaning of the YY-functional itself, 
it is sufficient for the calculation of the YY-functions, i.e., the critical values of YY-functional 
corresponding to the Bethe states. In this work we restrict our attention to the fc-vacuum 
state | v^™^ ) g Tik- In Section[3] it is shown that the corresponding YY-function can be 
identified with the on-shell action for the MShG equation. Another purpose of Section[3] is 
to discuss technical tools for the calculation of the YY-function. We review the well-known 
approach [2TII22] which allows one to express partial derivatives of the YY-function in terms of 
a solution to the non-linear integral equation from Ref. [2~2l . 

SectionH] is devoted to the so-called minisuperspace approximation (in the stringy terminol- 
ogy). The approximation implies the £ — > limit such that the soliton mass M is kept fixed. In 
this case the sine-Gordon QFT reduces to the quantum mechanical problem of particle in the 
cosine potential^] At the minisuperspace limit the world sheet of the MShG equation collapses 
into a single ray (see Fig. CP and the solution r\ at the segment (A, B) is expressed in terms of 
a solution of the Painleve III equation ( 11. 7p . subject of the boundary conditions 

f 2nog(*) + 0(l) as t^O 

1 — log(r - t) + 0(1) as t -> r 

Real solutions of the Painleve III equation which are regular at the open segment t G (0,r), 
and satisfy the boundary conditions (11.141) . form a family which is parameterized by r > and 



1 Note that in the conventional classical limit, the mass of the lightest particle in (jl.ljl is kept fixed while 
M -> oo. 
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(a) 



(b) 



Figure 2: The w-image of the cutted cone under the map (12. lj) . The points on the cone and 
their images are denoted by the same symbols. The segment AB is identified with AB, and 
the boundary line from B to infinity is identified with the line from B to infinity. The point O 
is an origin of the w-plane. 

— | < / < 2- By taking the minisuperspace limit of the generalized FSZ relation, we solve the 
connection problem for the local expansions of the solution U (t) (I1.14p at t — and t = r. The 
results obtained in SectionH] provide an interesting link between the Painleve III and Mathieu 
equations. 

2 On-shell action for the ShG equation 
2.1 From MShG to ShG 

In practical calculations it is convenient to trade the world sheet variable z in the MShG 
equation (ll.lOp to 

i7r(q + l) f , 

w = e 2 a I d z a/p(z) , (2.1) 

and similarly for w. The branch of the multivalued function (12 .ip can be chosen to provide 
the map of the cone with the cut along the ray (AB) visualized in FigfS^i to the domain of the 
u;-complex plane in Figf2b (see Ref. [TD] for details). This conformal map brings the MShG 
equation to the conventional Sinh-Gordon (ShG) form 

dvdvij - e 2f > + e~ 2f > = (2.2) 

for i) = 77 — I log (p(z)p(z) ), which vanishes at infinity 

lim fj = , (2.3) 

\w\— >oo 

becoming singular at the apex 

77 = 21 \og\w - w A \ +0(1) as \w — w A \ — >■ (2.4) 
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and at the point B ~ B 

f] — —— log \w — wb\ + 0(1) as \w — wb\ — > . (2.5) 
Unlike the apex singularity, the asymptotic (12. 5p is an artifact of the coordinate transformation 



2.2 Action functional 

To generalize relations (II .5p . (I1.9P we need an extra ingredient - the "on-shell" action for the 
ShG equation (12. 2p . It can be defined through the following limiting procedure. Start with the 
domain D depicted in Figj2b of the complex w-plane. Cut out the small sectors of radius e 
around the point A, B and B to obtain the domain D 6 shown in Figj3j Define the regularized 




Figure 3: The integration domain D e for the regularized action (12.61) . 



action functional 



A[fj[ 



lim 

e-S-0 



6ne 



— — — ( OwV<3wV + 4 smh (77) 1 H / d£ rj log(e) 

D e 2vri vre a 



d£ 7) 



1 



67T6 



di ^ ~ 12 log ^ e ^ 



(2.6) 



The first term is the "cutoff" version of the naive action for the ShG equation (12. 2p . The 
additional terms involves integrals over three arcs Ca, Cb and Cg and field-independent coun- 
terterms which provide an existence of the limit. Then the ShG equation supplemented by 
asymptotic behaviors near the singularities ( 12. 4ft . (12. 5 j) and at large u> ( 12. 3ft constitute a suffi- 
cient condition for an extremum of the functional (12.61): 



<u = . 



(2.7) 
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Finally we define the on-shell action A* as the value of A[fj] calculated on the solution rj 

For the variation (12.71) the world sheet geometry, as well as the parameter / controlling the 
behavior of the solution at the apex, is assumed to be fixed. Varying the on-shell action with 
respect to the parameter /, it is observed that 

BA*\ 1 

' -VA, (2.8) 



dl J r,a a 

where the constant t)a can be thought of as a regularized value of the solution fj at the apex 

fj A = lim (f)(w, w) — 21 log \w — Wa\ ) ■ (2.9) 

\w—wa\-^0 

It should be stressed that unlike I, which is the "input" parameter applied with the problem, 
the value of the constant tja is not prescribed in advance but determined through the solution, 
i.e. it is rather part of the "output". 

Let us consider now the infinitesimal variations of the world-sheet geometry. The corre- 
sponding 8 A do not vanish on-shell and can be expressed through the on-shell values of the 
stress-energy tensor. Under the infinitesimal dilation — = — = A <C 1, 



Sr 

8 r A* = — 
r 



f dw A dw n /I 2 1 
km / : 9 - - + — 

e^O J D 7T1 Va 12 



(2.10) 



where 

= 4sinh 2 (r/) (2.11) 

is a trace of the stress-energy tensor for the classical ShG equation. The other two non- vanishing 
components of T^ u are given by 

T = {8 W 1 1 )\ f=(^r)) 2 . (2.12) 

By virtue of the ShG equation, they satisfy the continuity equations 

d^T = d w e, d w f = d & e, (2.13) 

and, hence, they can be expressed in terms of a single scalar potential 

T = «9 2 $, T = a 2 $, Q = d w d^. (2.14) 
Combining the last formula with f )2.10p . one has 

/m-x =l . m r *»**> lKAfi _tP + i\ (215) 

V dr J a ,i e -^°JD. 71-1 va 12/ 

The 2-fold integral here can be reduced to the linear integral over the boundary of D e . The 
linear integrals over the arcs Ca, Cb and Cg from Fig J3] cancel out the term in the brackets in 
(I2.15p . This follows from the asymptotic formulas 

-2l 2 \og\w - w A \ +0(1) , \w-w A \^0 (2.16) 
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and 

Hw ^ )= l x (tog\w-w B \ + 0(l) \w-w B \^0 
18 I log \w — Wg\ + 0(1) , \w — w B \ — >0 

which are consequences of Eqs. (12.12p . (12.14ft and (12.41) . (12. 5ft . To proceed further, we need to 
use some properties of the potential $ discussed in AppendixEl Namely, for \w\ > \w B \ 

. hrja+l) i7r(a+l) . 

§(we <* , we ° ) = $(w, w) (2.18) 

and 

lim $(w, w) = . (2.19) 

Eq. (l2.18p implies that the half-infinite boundary rays (B, oo) and (5, oo) from Figj3] do not 
contribute into the integral (I2.15p . Now, taking into account Eq. (l2. 19j) . it is straightforward to 
show that 

r (^L=-i siii (^)( 3 > +5 ')- (22 °) 

where notations from Ref. [10] are adopted, 

sin(£)3i = i / (dwT + dwQ) (2.21) 

sin(£)3i = \ e~ iJ ^ [ (dwf + dwQ). 

Jc 

The integration contour C is visualized in Fig.HJ Due to the continuity equations, 3i and 3i are 
integrals of motion, i.e., they do not change under continuous deformations of the integration 
contour. 

Finally, let us consider the variation of the on-shell ShG action under an infinitesimal change 
of the apex angle -. In this case, using the simple electrostatic analogy, one can express S a A* 
through the torque applied to the boundary dD e 



5 a A* = S(-) lim [ - / — e^x^rf T va + -fj A + - log(e 

\OL/ c->0 |_ Jq d 7T 71 71 



2 

(2.22) 



where x 1 = dte(w — wa) and x 2 = Qm(w — wa) are real coordinates on D e , n u is a unit 
external normal to the boundary dD e and T^ u — — ~ d^f/d^r) + 5^ [ | (d a fj) 2 + 2 sinh 2 ^) ] . The 
integration contour dD t contains two components dD+ and dD~, related by reflection on the 
axis x 2 = 0. Since each component contributes equally, we replace the integral in (I2.22p . by 
2 J dD + and then evaluate it using the identity 

4 f dt e^x^rfT^ = [ dx"d^- [ d£ d^afQ) (2.23) 
Jc Jc Jc 
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t 



r/B 





Figure 4: The integration contour C in (I2.2ip . The contour on the cone and its w- image are 
denoted by the same symbol. 



and Eq. fl2~T9|) . This yields 



a 



OA* 



\<$>A-l f]A 



(2.24) 



da ) r,i 

where t)a is given by Eq. (l2.9p and $,4 stands for another "output" constant determined through 
the solution of the ShG equation - the regularized value of the potential at the apex 



$a = lim ( w) + 2l 2 log \w - w A \ ) • 

\w — Wa\— >0 

2.3 Generalized FSZ relations 

The compatibility of the derived above equations (12. 8ft . ( 12.2011 and ( 12.24jl implies 



(2.25) 



a 



03 



dl 
\daJ r,i 

where we introduce the notation 



- r (^) 

V or / a,l 

or ' ' - ' ■ " 



(2.26) 



dr J a,i 



5 = ^ sin (^)( 3l+51 )- 

According to Ref. [10] this constant is related to the sine-Gordon fc-vacuum energy 

R 



3 



7T 



Ek — Coo R ) , 



(2.27) 



(2.28) 



provided I = 2\k\ — |, a = £ 1 , and stands for the specific energy of the system with the 
infinitely large space size [23] : 

M 2 



lim 



Ek 



— tan — 

4 V 2 



(2.29) 



Thus, for < k < 5, relations (I2.26j) are recast into the form 

7T V <9£ 7r,fc ~ 8cosh 2 (^) 2 r V 9r ) a ,i V \ dr )a,l ' 

Formulas (I2.30p generalize the FSZ relations ( 11. 5B and ( II. 9p to arbitrary values of the sine- 
Gordon coupling constant and the quasi-momentum. Indeed, as £ = 2, the apex angle of the 
cone in FigfJ^ becomes 2ti, whereas k — j corresponds to / = 0, i.e., the solution of the (M)ShG 
equation remains finite at the tip A. In this special case, r](w,w) is expressed in terms of the 
Painleve III transcendent (jl.7p . (II. 8ft : 

f)(w, w) = U(A\w -w B \) . (2.31) 

Since \wa — wb\ — r/4 (see FigfSb), the value fj at the apex is given by 

Va = U(r) , (2.32) 

whereas, as it follows from the general relations d w du,& = 4 sinh 2 (r/) and (12.191) . 

I 

dt t sinh 2 U(t) . (2.33) 



d$A 
dr 



2.4 Normalized on-shell action 

Although the on-shell action A* disappears from the generalized FSZ relations, it is a main 
player in the derivation of (12. 30 p . Let us discuss some of its properties. 

The R.H.S. of (I2.28P exponentially decays at r — > 00 (see e.g. [21]). This enables us to 
represent the on-shell action in the form 

r°° dr 

A* = At + - $ , (2.34) 
J r r 

where the integration constant stands for lim^oo A*. The calculations outlined in AppendixlBl 
yield its explicit form 

A* x = log (3&2-§) + (3f + l) log (A G 2-1) (2.35) 

where Aq = 1.28243 ... is Glaisher's constant and we use the sine-Gordon variables £ = a^ 1 
and k = (21 + l)/4 > 0. 

The second term in Eq. 02.34]) is of primary interest thus we introduce the special notation 

POO J D 

2)= / — (E k - eoo R) . (2.36) 
Jr tt 
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Evidently it is the on-shell ShG action normalized by the condition 

lim 2) = . (2.37) 

1 — >oo 

Then Eqs.(ESD, (E2D are replaced by 



oo 



where we still assume that < k < 

Using the relation (see the conformal perturbation theory expansion (14. ip bellow) 

lim RE k = -~ c eff , (2.39) 
R^O 6 

where 

Ceff = 1 - y^y (2.40) 

is the "effective" central charge, one can represent 2) in the form which is appropriate for the 
study of the R — > limit, 

9 = | « log(Mfl) + ?)„ - ^ tan (|) _ f ( B4 + ^ ) . (j.4!) 

Here 2)o is some ^-independent constant. To calculate this constant explicitly one should write 
2J as A* — A^, express the functional (12.61) in terms of the original variables of the MShG 
equation (11.101) . and then analyze the limit of a small r. The straightforward calculations yield 

2j = T2 log(4^(l + 0- W )-|c eff logfo) (2.42) 



oo 



dx ( sinh(x) cosh(4£fcx) 1 c e ff _ 2x 



o 



2x sinh(£x) sinh (x(l + f )) 2 £( x + ^ 2 6 
where is given by Eq.f ll.12] ). 

3 YY-function in the sine-Gordon model 

In this section we identify 2) (I2.36P with the YY-function and briefly review the approach to 
numerical calculation of its partial derivatives. 

3.1 YY-function for the inhomogeneous 6-vertex model 

The sine-Gordon model admits an integrable lattice regularization based on the conventional 
i?-matrix of the six- vertex model (see Fig. [5]). Here I shall recall some basic facts concerning 
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R kl (2iQ) R?/(0) R%(-2$» 



Figure 5: Partition function = Tr[q k ^ a ? r^] of the inhomogeneous 6- vertex model on 

an infinite cylinder. Here r is the monodromy matrix along the infinite direction and q = e x +« . 
Rij(X) are conventional Boltzmann weights for the 6- vertex model satisfying the Yang-Baxter 
equation. 



the lattice BA equations which are relevant for the purposes of this work. All the details can 
be found in Refs. [lfl[25]. 

The energy-momentum spectrum in the lattice theory can be calculated by means of the 
algebraic BA, or Quantum Inverse Scattering Method: BA state is identified by an unordered 
finite set of distinct, generally complex numbers 9j which satisfy BA equations 



sifij + e + f ) a { 



e + f) 



?i + e-f) 



N 



4in£k 

-e !+« 



n 



n + m) 



0,, 



m) 



(3.1) 



where 



s(x) = sinh ^ 



x 



(3.2) 



and iV stands for one-half of the number of sites along the compactified direction in Fig.0 
The parameter O controls the world-sheet inhomogeneity of the Boltzmann weights, whereas 
k in ( 13. ip is proportional to the twist angle for the quasiperiodic boundary conditions imposed 
along the compactified direction. Then the energy E^ N ' and momentum of the BA state 
can be extracted from the formulas 



exp 



E (N) ± p(N) 

2N 



n 



s(f + e ± 6j) 



(3.3) 



For the vacuum state all the BA roots are real and their number coincides with N, which is 
assumed to be even bellow. Following Yang and Yang |16j . the BA equations in this case can 
be bring to the form of the extremum condition 



o 



N 
2 



N 
2 



1) 



(3.4) 
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for the YY-functional defined by the formulas: 



with 



Y w) = 2 ( m) - ) + £ ^ - <w < 3 - 5 > 



n^-^r ^'-lrf' e" (3-6) 



7T 7_oo sin h(Mi±a 



(3,) 



and 



7-oc S inh(^i±^ 
Here and bellow the symbol f stands for a principal value integral defined as the half-sum 

1 i r+oo— iO p+oo+iO \ 

2 V J-oo -iO "J " J-oo+iO /' 

Eqs. (13.41) can be interpreted as an equilibrium condition for the system of N one-dimensional 
"electrons" in a presence of confining and linear external potentials. For large separations, 
0»1, 

U ^ = ~Y^ \0\+O(e-$&) , (3.8) 

therefore the 2-body potential is essentially a ID repulsive Coulomb potential slightly modified 
at short distances. Since 

Nf N£ , 2|fl±e| , 

^» = 2(rTi)i 9 - e i + 2(rTi)i (,+e i +0 ( e "^)' < M > 

V{6) can be interpreted as a potential produced by two heavy positive charges + 2 n+ti pl ace d 
at ±0. We shall always assume that the external linear potential in ( 13.51) is sufficiently week 
and the inequality 

-\<k<\ (3.10) 

is fulfilled. For 

0<£<1, (3.11) 

the Hessian of the system (13.41) . , is positive definite therefore we will focus primarily on 

this case. 

As the physical intuition suggests, the YY-functional (13. 5 j) has a stable minimum at some 
real distribution of the BA roots 



e {N l < ew +1 <...< of_ 3 < ef_, . (3.i2) 
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Figure 6: from Eq.f l3.15|) for A" = 400 and = 2, A; = 0, £ = |. The solid line represents 
the continuous density A3. 16|) . 

The main subject of our interest is the YY-function, i.e., a critical value of 7^' calculated 
at this minimum. With some abuse of notations we will denote it by the same symbol as the 
YY-functional, 



(3.13) 



Using the YY-function, the ground state energy (13. 3p can be written as 




(3.14) 



whereas the momentum associated with the ground state is of course zero. 
At large A" and finite 0, the distribution of the BA roots 





(3.15) 



is well approximated by the continuous density (see Figj6]) 



P(0) 



l r l l 



(3.16) 



2tt Lcosh(#-0) cosh(0 + 6) 



Therefore the following limit 



does exist 



lim N- 2 Y^ N \Q) = yoo (e) 

N— >oo 



(3.17) 



and it is a simple exercise to 



show that 




(3.18) 
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For large 9, Eq.f l3.18p yields 



*»(e) = |g + + 1 ^ tan (f ) + (.-») . (3.19) 

The constant term here coincides with half of the value of (13.181) taken at = 0. This is not 
an accidental relation. Indeed, as |0| — » +oo, all the BA roots split into two clusters centered 
at ±0. The systems of BA equations for each cluster are completely separated in this limit 
and reduce to the original form (13. ip with = and N is replaced by iV — > N/2. Hence for 
any even N 

y(*)(e) = i^l |0|+2F W2) (O) + o(l) as -> ±oo , (3.20) 

where the first term describes monopole-monopole interaction of the "electron" clusters while 
the second one represents their intrinsic potential energy^] Combining (13.14j) with (I3.20p one 
also has 

S w (6) - = o(l) as 0^+oo. (3.21) 

It should be emphasized that asymptotic formulas (I3.20|) and (13.211) do not assume the large- N 
limit and can be applied for any finite N. 

3.2 Scaling limit 

The sine-Gordon QFT (II. ip manifests itself in the scaling limit when both N, — > +oo while 
the scaling parameter 

r = ANe- @ (3.22) 

is kept fixed (RG-invariant). In this case the L.H.S. of ( 13 . 2 1 [) does not vanish, but has a simple 
relation to the fc-vacuum energy 



lim ( - 1^ "I = ™* + £* , (3.23) 
- • 1 1 + f J 2tt 12 ' K ' 



r — fixed 



where the effective center charge c e g is given by Eq. (l2.40p . 

In order to study the scaling behavior of the YY-function, it makes sense to consider only 
the part of the "electron-ion" potential energy corresponding to the mutual interaction of the 
clusters, 

Y£P(Q) = Y (7V) (0) - |^ |0| - 2 Y^ 2 \0) , (3.24) 



2 In a view of the mechanical analogy, it would be natural to include an additional term — |0| into 
the R.H.S. of definition (|3.5[) . This term represents the ion-ion potential energy and does not affect on the 



equilibrium conditions (|3.4[) . 
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which vanishes as — » +00 for any fixed N. Taking into account Eqs.f l3.14j) . (13. 23 p . we get 



r -fixed JU 

or, equivalently using 

lim ^(O) = 2) - i c efr log(r) + £ tan (Zl) - 2) - (3.26) 



N,e- 

r — fixed 



The last formula allows one to identify the YY-function in the sine-Gordon QFT with the 
normalized on-shell action 2). 

The following comment is in order here. Our analysis is based on the existence of solutions 
(I3.12p of the vacuum BA equations. It can be directly applied to the case < £ < 1 only. At 
£ = 1, the sine-Gordon model is equivalent to the free massive Dirac fermions theory and a 
closed form of the YY-function can be easily derived from definition (I2.36H : 

/°° dr r 1 

— r sinh(r) log (l + e - rcosll M+ 27rife )(l + e -rcosh(r)-2 7r ife^ _ ^_ 2 ^ 
-00 2tt L 

It is expected that for £ > 1 the YY-function is uniquely defined through the analytic contin- 
uation from the segment £ G (0, 1). 



3.3 BA roots at the large-iV limit 

Properties of the BA roots at the scaling limit were discussed in Ref . [25J . The roots accumulate 
at 9 = ±0. However, at the center region (see Figj6]) and at the tails of the distribution the 
roots remain isolated and their behavior can be described as follows (see TablesQ] and [2] for 
illustration) : 

• There exist limits 

6 3 = lim 9f (j = 0, ±1 ±2...) (3.28) 

iv,|0|— too J 
r, j — fixed 

and, for an arbitrary > 0, 

ri +) = lim (e {N) N +\og(N) + Q) (3.29) 

n — fixed 

7-H = lim log(iV)-e) (n = 0,l...). 

m — fixed 

3 Notice that (0) can be interpreted as the YY-function for the spin-i Hciscnberg chain. As N — > 00 

Y (N \0) = ^(0) N 2 + I c cff log +2) + , 
where 2)o and z/oo(0) are given by Eqs. (|2.42p and (|3.18[) . respectively. 
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n 


N = 100 


N = 200 


N = 400 


N = 800 


iV = 1600 


a 
U 


i a a o /i o 


1 a /i o /i o 


1 A/I Q /I A 

l.U4o4U 


l.U4o4U 


1 A/1 Q A A 

l.U4o4U 


1 


3.01807 


3.01640 


3.01598 


3.01588 


3.01585 


2 


5.01975 


5.01202 


5.01009 


5.00960 


5.00948 


3 


7.03507 


7.01378 


7.00848 


7.00716 


7.00683 


4 


9.06565 


9.02023 


9.00896 


9.00615 


9.00545 


5 


11.1150 


11.0317 


11.0111 


11.0059 


11.0046 


6 


13.1873 


13.0489 


13.0149 


13.0064 


13.0043 


7 


15.2870 


15.0729 


15.0204 


15.0074 


15.0042 


8 


17.4186 


17.1044 


17.0280 


17.0090 


17.0043 


9 


19.5874 


19.1447 


19.0378 


19.0112 


19.0046 



Table 1: ^ exp (9 { n N >) for r = 1, ^ = | and k = 



n 


N= 100 


iV = 200 


N = 400 


N = 800 


N = 1600 





1.02837 


1.02831 


1.02830 


1.02830 


1.028299 


1 


3.01276 


3.01111 


3.01069 


3.01059 


3.01056 


2 


5.01652 


5.00881 


5.00689 


5.00641 


5.00628 


3 


7.03273 


7.01148 


7.00619 


7.00487 


7.00454 


4 


9.06381 


9.01843 


9.00718 


9.00436 


9.00366 


5 


11.1135 


11.0302 


11.0096 


11.0045 


11.0032 


6 


13.1860 


13.0477 


13.0136 


13.0051 


13.0030 


7 


15.2858 


15.0718 


15.0194 


15.0063 


15.0031 


8 


17.4176 


17.1035 


17.0271 


17.0081 


17.0033 


9 


19.5864 


19.1438 


19.0369 


19.0104 


19.0038 



Table 2: ^ exp + ©) for r = 1, f = § and k = 
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The limiting values of the roots possess the following n — > +00 asymptotic behavior 

e dn = — (2n + l + 2k) +0(n- 1 ) (3.30) 

27T 

e -8-n-i = _ (2n + l -2k) +0(n- 1 ) 



and 



exp(±^ ±) ) = | (2n + l±2k) + 0(r^ 1 ) . (3.31) 



To probe the infinite sequences {0j}°Z_ oo and {ri }™ =0 it is useful to consider certain generating 
functions encoding their properties. Let (+(ou) and C-( w ) De functions defined as the analytic 
continuation of convergent series 



00 



= ^ e -^„ Sm(u) < -1 , 

ra=0 
00 

£_( w ) = ^ e -i^-i-« 3m(w) > 1 . (3.32) 



n=0 



As follows from the asymptotic formulas O3.30p . (+(uj) (£_(u;)) is analytic in the half plane 
^sm(u) < 1 (3m(u) > —1) except at a simple pole at uj = — i (u = i) with the residue — f- 
Also note that C±(0) = =F/c. Therefore 

C(w) = C+H + C-(w) , (3.33) 

is an analytic function in the strip |Q : m(a;)| < 1 such that 

C(0) = . (3.34) 

In the limit r — > 

C(«) = (£f (CrV) + o(i)) + -,W (Cif(-o;) + o(i)) , (3.35) 

where (±o;) are the zeta functions for the sequences {ri }^L 1 )3. 29ft . i.e., 

00 

Cif(±^) = E ex P(- i ^ ±) ) • ( 3 - 36 ) 

n=0 

The function ^ cft ' ) (a;) was introduced (in a different overall normalization) and studied in Ref. 
[26J. Its exhaustive description was found later in Ref. [27] (see also related works [28] and [29]). 
where it was shown that Cfc ( w ) coincides with the zeta function of the Schrodinger operator 

- d 2 x + x 2a + ^4r^ . (3.37) 

ar 
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More precisely, for < k < |, the sequences {r,i are simply related to the spectral sets 

[£^}%L of this differential operator: 



exp 



(3.38) 



with a = I = 2k — | and given by (11.12p . 

With the above properties of the BA roots it is not difficult to analyze the large iV-limit of 
the relation 



\ dk )n,b& 1 + f ^ ■ 



(N) 



(3.39) 



which is derived by differentiating (13. 5 j) with the use of BA equations (13. 4p . The scaling analog 
of (13T39P reads 



\dk)r£ 1 + f 



4i( C'(0) , 



(3.40) 



where prime stands for the derivative with respect to u. As follows from (I3.35P and Q (0) = 

-k, 

(f = -^f ,og (D - ni (^'» + a "^ ,(0) ) + as - • < 3 ' 41 ' 

The subleading term of this asymptotic is expressed in terms of the determinant of the differ- 
ential operator (I3.37P and can be calculated explicitly: 



( Op / r,£ 



2p 



^ r(i + |)r(i-^; 



r(i-f)r(i + I ^; 



+ o(l) as r -> 0, (3.42) 



where we substitute for the equivalent parameter p = 2£k. Of course, Eq. (l3.42p can be 
alternatively obtained by means of relations (I2.4ip . (I2.42p . Note the expression in the bracket 
[ • • • 1 coincides with Liouville reflection amplitude (analytically continued to the domain — 1 < 
£ < 0) introduced in the work 



3.4 Calculation of partial derivatives of the YY-function 

For < £ < 1 the so-called Q-function can be defined through the convergent product 



Q{6) = £ e^t []4eT s (0 B -(l) s (9-e_ fl _ 1 ) 



(3.43) 



n=0 



where the abbreviation s(6) (13. 2p is applied. The ^-independent factor £ can be chosen at will. 
In what follows it is assumed that 



e = 2^ H e ^\m ( 21 °g(2^TT))+^-^-i) 



n=0 



(3.44) 
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In the scaling limit, BA equations (I3.4p boil down to 

e (^) = vr(2j + l) (7 = 0, ±1, ±2... 

where 



and the branch of the log is fixed by the condition 



(3.45) 



(3.46) 



re 



e{9) = — - 2«k + 



7T 



for Ke(6>) ->■ +oo and |9fm(0)| < - . 



(3.47) 



Using the analytic properties of (±(uj) fl3.32[) . definitions (I3.43P and (I3.46P can be trans- 
formed into the integral representations 



logQ(# + i^) = mk + 



r cosh(#) 1 
2cos(f) ~ 2 



du 



CM 



w sin h(^i±i); 



(3.48) 



and 



e(0) = -2vrA; + r sinh(0) - i 



du sinh( 7raJ ^ ^ 



C(w)e 



respectively. On the other hand, the BA equations ( I3.45P imply (see Ref. 
Ref. [21] for the original derivation) 



CM 



21 > 



d0 



log (1 + 



-ie(0-iO) 



(3.49) 

and related 
(3.50) 



cosh(^f) sinh(^) J-oo 2tt 
Note that at the free-fermion point (£ = 1) e(8) = r sinh(0) — 2irk, therefore Eq. fl3.50p gives 
d0 



e 2 log ( 1 



-r cosh(0) +27rifc 



) -e^ \og(l + e- rcosh ^ ik ) . (3.51) 



In general, Eqs. fl3.49p and (I3.50p are combined into a single integral equation on e(8). Once 
the numerical data for e(9) are available, ((uj) can be computed by means of (13.501) . 

Eq. fl3.50p shows that ({uj) is a meromorphic function with simple poles located at uj = 
±i (2n+l) and uj = ±j (n+1) (n = 0, 1 . . .). The residue values are the k- vacuum eigenvalues of 
local and nonlocal integral of motions in the quantum sine-Gordon model [T01I26"] . In particular, 
at the boundary of the strip of analyticity |Sm(o;)| < 1, ((uj) has simple poles 



CM = T-cot (f 
r \ 2 

where $ is given by Eq. (I2.28|) . 



uj ± i 



+ 0(1) as uj ->■ ±i 



(3.52) 



20 



This way, the problem of numerical calculation of k- and r-partial derivatives of the YY- 
function can be solved by means of relations 

= tan(^) Um.(l=Fiw)C(w) , 



V dr ) i,k 



\dkJr,a 1 + f 



C'(0) . (3.53) 



The calculation of ^-derivative is found out to be a more delicate problem. Rather naive 
manipulations with the lattice YY-functional (13. 5p suggest that 

/<92J\ 1 r° du sinh(Trw) ... . fnrA . 

In AppendixO we present some evidences in support of this relation. Unfortunately it still 
lacks a rigorous proof. 



4 Minisuperspace limit 

4.1 Minisuperspace limit of the YY- function 

The small- R expansion of the /c-vacuum energy in the quantum sine-Gordon model was argued 
in Ref. 0, 

RE k 1 Ak 2 S ^ 4n 

7T 1 + C Z — ' 

n=l 

The first coefficient f\ has a relatively simple explicit form 

/ g(l-2fc) x / g(l+2fc) x 

/i =4 T(l+ f ( ^ ^ J , (4.2) 

where 7(2) = T(x)/T(l — x) and (11.121) . Let us consider the £ — >• limit of (14. lft in which 
the parameters r and are kept fixed. One finds 

— = 4+^ + «(0 (4-3) 

with 

a = " 2 + 2(^1) + ° (g4) - (44) 



Here we have denoted 



g=(^)\ z/ = 2A;. (4.5) 
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In this special (minisuperspace) limit the sine-Gordon QFT reduces to the quantum mechanical 
problem of particle in the cosine potential whose energy coincides with a from Eq. fl4.3p FI Note 
that ( 14. 5 p are conventional notations in the theory of Mathieu equation [32] . For given is, a is 
determined by the Whittaker equation 



• 2 (™\ 



sm ~ lyj = A? ( a ) sin2 ( 7rx >" 



where A 9 (a) is Hill's determinant 



A g (a) = det 



V 



• • ■ 7-2 1 


72 




7o 


1 


7o 




72 


1 72 • • • 



\ 



/ 



72r 



An 2 



(4.6) 



(4.7) 



The solution of Eq. (14. 6J is a multivalued function, but the condition ( 14. 4 p specifies the proper 
branch unambiguously!^ To simplify formulas, we will below treat a as a function of the variables 
r and is, i.e., 



a = a(r, is) . 



(4- 



Now, using Eqs.f l2.4TT) and ( 14. 3p it is straightforward to obtain the limiting behavior of the 
YY-function: 



where 



y(r, is) 



with 



00 dx 
x 



sinh (isx) 



v 



iq ju y xsmh{x 
Note that the above form of y(r, is) can be alternatively rewritten a: 



y(r, is) = log (2« A%) + - log(r) - T - 



°°dt ( . , t 2 t 1 

— a(t,u) H h - 

t V v ' ; 8 2 4 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



4 The minisuperspace approximation for the closely related quantum sinh-Gordon model was discussed in 
Ref. [31]. 

5 It is implemented in the Mathematica as MathieuCharacteristicA[ v, q] with — 1 < v < 1. 
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For the large-r expansion of a see, e.g., formulas 28.8.1, 28.8.2 in Ref. [34] . 



22 



4.2 Minisuperspace limit of the Q-function 

The minisuperspace approximation of the Q-function at r = was argued in Appendix B 
of Ref. |26j. The analysis was based on general properties of the Q-function and only minor 
adjustments are needed to extend it to the case r > 0. 

The Q- function is a quasiperiodic solution (see Eq.f l3.43p ) 

Q(0 + i7r(l+O) =e 2ink Q(0) (4.13) 
of Baxter's equation (see, e.g., Ref. [TO] ) 

Ti{9)Q{9) = Q{9 + mi) + Q(9-m£) , (4.14) 

where Ti (9) stands for the k- vacuum eigenvalue of the transfer matrix. If the overall normal- 
ization factor in (13.431) is chosen as in Eq. (13.441) . then the Q-function also obeys the so-called 
quantum Wronskian relation 

Q(0 + if, k)Q(9-'^,-k) -Q(9- '^, k)Q(9 + ^,-k) = 2i sin(27rfc) , (4.15) 

where we explicitly indicate the dependence on the quasi-momentum. In the minisuperspace 
limit 

Ti (9) = 2 + «) 2 w(9) + o(0 , (4.16) 

whereas Baxter's equation reduces to the second order differential equation. For the conformal 
case (i.e., for r = 0) discussed in [26], w(9) = e 2e — {2k) 2 . Similarly, in the case of finite r one 
can show that 

,2 

w (9) = L CO sh(2^) - A , (4.17) 
8 

with some ^-independent constant A such that 

lim A = (2k) 2 . (4.18) 

r- >0 

Thus the minisuperspace limit of the Q-function can be described as follows. Let F u (z) be 
Floquet's solution 

F u {z + m) = e 1 ^ F u {z) , F- V (z) = F u {-z) (4.19) 

of the modified Mathieu equation 

d 2 F ( r 2 \ 
- + (A - — cosh(2^) J F = (4.20) 

normalized by the condition 

W[F v ,F-„] = -2 sin(7rz/) , (4.21) 
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where W[f, g] stands for the Wronskian fg' — gf. Then, as follows from (I4.13p - fl4.17p . 

lim Q(6, k) = F v {6) with v = 2k . (4.22) 



, r. k — fixed 



For given u, the constant A in (I4.20p is determined by the quasiperiodicity condition (I4.19p . 
which implies the Whittaker equation (14. 6p with a replaced by A. The extra condition ( 14 . 1 8 j) 
enables us to chose the branch of solution of (14. 6 p unambiguously. Thus we conclude that 

A = a(r, v) , (4.23) 

where a is the same function as in Eqs. fl4.3p . (14. 8p . It is useful to note that Eq. fl4.2ip is equivalent 
to the following normalization condition^ 



dy F v (iy) F v (-iy) = 27rsin(7rz/) 



/da 



(4.24) 



Eq.f l4.22p dictates that the BA roots {9j}°? = _ 00 (I3.28P in the minisuper space limit turn out 
to be the zeros of the Mathieu functional 

lim Oj = zj : F v (zj) = (j = 0, ±1, ±2 . . . ) . (4.25) 

j. r, k — fixed 

Using the properties of F u (z), it is not difficult to derive the sum rule 

r°° dx 

S = Sin M / rp ( ,^JP(_ _»rN » ( 4 - 26 ) 



where S stands for the regularized sum of (—2Zj) 



S = 2 £ ( - *n - z-n-i ) ~ 2, log (^) (4.27) 



n=0 



{pfE is Euler's constant). Note that £ = — 2i lim r _> C'(0), and hence Eqs. (13.401) . (14. 9p imply 
that 



dt /da\ 
T \dv)t 



(4.28) 



7 For < v < 1, F v (z) = M v (ce„(i;z) — \se v (\z) ) , M v = se f^cr jo) > where ce v (x) and se v (x) are returned 
by the Mathematica functions MathieuC[a, q, x] and MathieuS[a, q, x], respectively. Their x-derivatives are 
implemented as MathieuCPrime[a, q, x] and MathieuSPrime[a, q, x]. Here a = a(r, v) and q = ( j) . 

8 The remaining terms in the large-n asymptotic formulas (|3.30[) diverge as £ — > 0. For this reason 

these formulas are not applicable at the minisuperspace limit. 
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4.3 Connection problem for the Painleve III equation 

We now turn to the ShG equation (12.21) at the minisuperspace limit. In this limiting situation 
the triangle ABB in Fig|2b shrinks to a segment while fj becomes a certain solution of the 
Painleve III equation (11.7)) 

lim 77(10, w) = U(A \w — w A \) (0 < \w — w A \ < r/4) , (4.29) 

r, I — fixed 

such that 

e 2U(t) = K 2 t 4»-2 + o( ^- 2) for t ^ q _ (4>30) 

Parameters v and k are related to / and tja (12.9)) as follows 

v = l + i (4.31) 
log k = —I log(4) + lim r\A ■ 



a — »-oo 
I, r — fixed 



Let us discuss the solution satisfying (14.30)) in the context of the general theory of Painleve 
III equation. For < v < 1 and real k, the asymptotic condition ( I4.30[) unambiguously 
specifies a two-parameter family of real solutions of the Painleve III equation. A systematic 
small-t expansion of ( 14 . 3 [) has the form of the double series [7] 



e 2E/(t) 



K 2 t 4u-2 + B m , n t 4 ^ m+Avn - 2 (4.32) 



m,n — 
Tra+n> 1 



whose coefficients are uniquely determined through the parameters v and k by a recursion 
relation which follows from the differential equation (II. 7p . 

The expansion ( 14 .32)) is expected to converge for sufficiently small t. Let t = r be the closest 
singularity to the origin. The differential equation (II. 7p possesses the Painleve property which 
states that, except at t — and t = 00, the only possible singularities of e 2U are the second 
order poles of the form e 2U ^ = 7^332 — r ^_ r ^ + C + o(l) with some constant C. Further terms 
in this Laurent expansion are expressed in terms of r and C. They can be easily generated 
through the differential equation (II. 7p . It will be convenient for us to substitute C for an 
equivalent parameter c such that 

nnt\ 4 4 13-16c 2(160-7). , x2 x 

(t — r) z r (t — r) 6r z 6r 6 

We will focus on the case when the closest pole to the origin is located at the positive real 
axis, i.e., r > 0. This requirement imposes certain constraints on admissible values of k and 
c. Within the admissible domains, each pair (u, k) or (r, c) can serve as an independent set of 
parameters for the two-parameter family of real solutions of the Painleve III equation which are 
regular at the segment t G (0, r) and characterized by the behaviour U — > (2v — 1) log(t) + 0(1) 
as t — > 0; However, it is more convenient to choose (r, v) with r > 0, < v < 1, as a basic set 
of independent parameters. 
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At this point, we turn to the problem of finding the functions K = k{t, v) and c = c(r, v), 
i.e., the connection problem for the local expansions (I4.32p and (14.331) . It is relatively easy to 
establish the relation 

^ . (4.34, 



du/r \ dr 

The proof is similar to our previous derivation of the generalized FSZ relations; One should 
consider the action functional 



S[U] = -lim 

L 1 4 e^O 



J dtt((j 2 + smh 2 (U) - (r ^ t)2 ) + 2(2z/-l)f/(e)-(2z/-l) 2 log(e) 



(4.35) 



where the dot stands for the t-derivative. For U(t), t e (0,r) satisfying the boundary condi- 
tions (14.301) . (14.331) . the functional S[U] is well defined and its variation vanishes provided U(t) 
satisfies the Painleve III equation and 5U(r) = 0. Let S* be the on-shell value of (14.35)) . One 
can show that 

r(^f-) = X -- r --c (4.36) 
V dr Ju 4 8 V ; 

dS*\ 

-7T- = log K , 

OV J r 

and the compatibility of these equations implies (I4.34j) . 

Now let us apply the results of the previous sections. At the minisuperspace limit the first 
generalized FSZ relation (12.301) yields the formula identical to ( 14. 34ft with c(r, v) replaced by 
a(r, v). Therefore, we conclude that c — a does not depend on v. It is straightforward to analyze 
r — > behaviour of the on-shell action: 

= ( \ ~ v 2 ) log(r) + v log (^) - \ log Q - ^ + 0(r 4 ) . (4.37) 

This asymptotic formula, combined with (I4.36j) . implies that c = z/ 2 + 0(r 4 ), and hence c — a = 
0(r 4 ) (see Eq. (l4.4p ). In Appendix|D] it is explained how to systematically recover the small-r 
expansion of c(r, v). The calculations yield the expansion 

2 . 1 /r\ 4 . 5^ 2 + 7 /r\8 



" + 2(^-1) liJ + 32(^-1)3(^-4) UJ (438) 
9 z/ 4 + 58 v 2 + 29 /r\i2 ^ lg 



+ 6 4( z/ 2_ 1 )5( z/ 2_ 4 )( z/2 _ 9 ) (,4 ) +°( r ■ 

which matches exactly the small-r expansion of a(r,u) (see Eq. (l4.4p and formula 20.3.15 in 
Ref. [32J). The formal derivation of the relation 

c = a(r, v) (4.39) 

can be obtained with the use of equations (1A.31) . (IA. 1 7[) from AppendixlAl combined with the 
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0.2 0.4 0.6 0.8 1.0 (b) 

Figure 7: The admissible parameter domain for the two-parameter family of solutions of the 
Painleve III equation: (a) in terms of (y, n) from (I4.30p . (b) in terms of (r, c) from (I4.33p . 



results from Section POln Eq. fl4.39p . together with 

dt 



K 



Iv r 2v exp 



>1— 2i> 



r(i-^) 



t 



exp 



2v 



t) 

.ov/t 
dt /da 
~t Vdu 



(4.40) 



leads to an explicit solution of the connection problem for local expansions (14 .321) and (I4.33p . 
The domain of applicability of Eqs. (14.391) . (I4.40H is given by the inequalities (see Fig.[7h) 



< v < 1 , 



K > 



)l-2u 



T(l 



V) 



It can be equivalently described in terms of the pair (r, c): 

r > , Oo( r ) < c < bi(r) 



(4.41) 



(4.42) 



where ao = a(r, 0) and hi = lim^i a(r, v) stand for the minimum and maximum of the first 
conduction band for the Mathieu equation, respectively (see Fig. [7b)- 

Finally, let us briefly discuss the minisuperspace limit for = | ( $ — fj). Contrary to 
r)(w, w), the potential $(w, w) does not have a finite limit for < \w — Wa\ < f/A. However, the 
divergent part of &~>(w, w) for < \w — wa\ < r/4 is somewhat trivial and can be resolved by 
by means of the decomposition (1A.13I) from Appendix|Al The non-trivial part is specified 
by the conditions (IA.14|) .( p\.15|) and remains finite as a — > oo. It is convenient to introduce 



W(t) = -/(/ + 1) log \w -w A \ + lim $ ( } 



w, w) 



a— too 
r, I — fixed 



(4.43) 



Note that Eq. (|4.39|) allows one to derive the following large-a asymptotic formula for the on-shell action 



3^ r 
2 4 a 



a 



ea ar 



))+o(a- 1 ) 
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where t = 4 \ w — wa\ < r and I = v — ~. As follows from Eqs. (1A.2j ). (1A.13j) . W(t) satisfies the 
linear inhomogeneous differential equation (assuming U is given) 



t dt\ dt J 

and the boundary condition 

W(t)= (\~v 2 ) log o(l) , as t->0. 



Therefore 



W{t) 



4 



- v 2 lo 



t\ t 2 



16 4 



i / dTTl0g ( 



-2f/(r) 



(4.44) 



(4.45) 



(4.46) 



The t-derivative of W(t) at t — r is given by 



rW(r) = 1 --^-u 2 + \ 



dtte~ 2U ^ 



and it is not difficult to show that 



W{r) 



dS* 
Or 



The last two relations combined with Eqs.f l4.36j) . (I4.39P imply 

i f dtte- 2U ^ = v 2 -a(r,u) 



(4.47) 



(4.48) 



(4.49) 



By transforming the integrals as in Ref. [35J, the value of W(t) at t — r can be expressed in 
terms of a, k and S*: 



W(r) = ~ log (~) - j + (y - + 2z/ 2 log(2) - a + v log(«) - 5* 



(4.50) 



which is equivalent to the following relations 



W{r) 



or 



W(r) = y {v)-log 



A 6 G e—*~ v 2^+3 



Vl-u) 



(4.52) 



+ 



dt / (da 



dt 



da 



T * "ST + " -a(*^) + T-7 



<9z/ 



t 2 1 



8 4 



where 3^0 (^) is given by Eq. (14. lip . Note that the last term in Eq. (14.521) vanishes as r — > oo 
while the remaining part reproduces the result quoted in Ref. [7]. 
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5 Concluding remark 



In this work we have described the link between the action functional for the classical ShG 
equation and the YY-function corresponding to the fc-vacuum states in the quantum sine- 
Gordon model. The natural question arises: Can this relation be generalized for the excited 
states? Nowadays the machinery of the Destri de Vega equation for the excited states are well 
developed [25j|36], so that the calculation of the YY-function for the excited states does not 
seem to be a particularly complicated problem. However, it remains unclear how to construct 
integrable classical equations associated with the excited states and perhaps more importantly, 
what all of this really means. 
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A Appendix: Basic properties of the potential 

Let us define the following linear combinations of the solution fj of the ShG equation (12. 2p and 
the potential $ fl2TT41 

$W = ± ($±r)) . (A.l) 

They satisfy the relations 

a*$(±> = I ( (d,^) 2 ±d 2 w fi), dl&V = \ ( {d^f ±dlfj) , flWofcW = ^ - 1 , (A.2) 

which can be considered as a closed system of nonlinear partial differential equations. Results 
of Refs. [7] and [10] imply that the desirable solution of (1A.2j) is expressed in terms of the 
Fredholm determinants 

$(±)( W! ^)=logdet(l±^^ . (A.3) 
The kernel of the integral operator reads 

cosh (-£-) 

g w ,M = + ex V (-2we d -2we- e ) , 

where Ti(9) is defined by the Q-function (I3.43P through Baxter's equation (14.141) with £ = a -1 . 
In ( 1A.3I) (w,w) play a role of complex parameters. Within the domain 

|argH|<^, (A.5) 
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(w , w) can be represented by convergent series 

(±1) 



oo 



$&(w,w) = V^^- $ n (w,w) , (A.6) 



n 

n=l 



/oo n 
n 
-oo -_i 



d9j g 



W,W \ V J J 



-oof = \ V 47r COsh(M±i) 

where it is implied that 9 n+ i = Q\. Some immediate consequences of (]A.6j) . 

-) = $W (w, u>) for \w\ > \wb\ (A. 7) 

and 

lim $ (± \w,w) = , (A.8) 

| -La | — >oo 

have been used in the derivation of basic Eq. fl2.20p . By deforming the integration contours in 
(IA.6I) . the applicable domain (IA.5[) can be slightly extended to the triangle OBB (see Figj2b). 
This leads to the fact that for $le(w — wb) > 

& ± \w,w) = & ± \w,w) . (A.9) 

Yet the expansions ( IA.6j) cannot be applied inside the triangle ABB, 

71 

< %le(w — wa) < $te{wB — wA , arg(w — wa) < tt - ■ (A. 10) 
Within this domain, fj(w,w) satisfies the symmetry relations 

v{ w a + (w ~ wa) e«" , wa + [w — wa) e - "" ) = fj(w, w) = fj(w, w) , (A. 11) 



whereas, as follows from ( 1A.2I) . 

d 2 w <J$(±) = dl = d w d a = , (A. 12) 

where 5$ (±) = $ (±) (w A + {w-w A ) e^, w a +(w-wa) e"^ ) -$ (±) (u>, w). From Eqs. (J23| , flAJj]) 
and the reflection symmetry relation (1A.9|) applied to the segment -B.B, we deduce that within 
the domain (IA.10I) $ ± can be written in the form 

$> {±) (w,w) = -1(1 =f 1) \og\w -w A \ + (w- w A )J+ (w - w a )J + ® { a ) + $ ± (w, u)) , (A.13) 

where $^ and J = J are some real constants, 

$^ ± - ) (w^ + (it; — wa) 6«" , + (w — iua) e~^ ) = $ ± (w, w) = $ ± (u), iu) (A. 14) 

and 

lim ¥ ±) (w,w)=0. (A.15) 

\w— wa\-*0 
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The constants in (1A.13I) are given by the linear combinations of fjA (12.91) and $a (12.2511 : 



(±) 



(A.16) 



As for J and J, the calculations performed in Section l2T2l imply tha 



D 



dw A dw i +2fl \ 1(1 =F 1) , 1 1 T . ( 7T \ / _ ^ N / » i «\ 

4(e ±2 "-l) = ^-^±- + --- sm(-) • (A.17) 
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a 



V2a/ 



On the other hand, the integrand here coincides with the Laplacian 4 9 u; 9 iD $^ ± - ) , and the L.H.S. 
of (1A.17|) can be expressed in terms of J = J. This yields the relations 



J = 3i 



(A.18) 



Finally, note that Eqs.( l2.5j ). flA.2[) yield the following asymptotic behaviours in the vicinity 
of B ~ B (see Fig. [2b) 



$ (+) ->■ -— log |w - w B \ + 0(1) 
3d 



$ ( ) -» — log |w 
36 61 



u;j|+0(l). (A.19) 



B Appendix: Calculation of = lim r ^ 00 .4 



For large r the dominating contributions to the on-shell action (I2.6P come from the vicinities of 
the points A and B ~ B (see Figf2b). Near these points, fj can be approximated by U(4\w — 
w a\)\ u =i + i an d U(A\w — wb\)\ u= ^ respectively, where U are the Painleve III transcendents, i.e., 
regular at t > solutions of Eq. (II. 7p satisfying the boundary conditions 



U ->■ 



(2z/ - 1) log (|) + l/oM + o(l) , as t-> 



, as t — )• oo 

This observation implies that the limiting value of the on-shell action ( 12. 6 j) is given by 



where <S^(^) is the on-shell value of the action functional for the Painleve III equation, 



(B.l) 



Soo[t/] = - Jim 



J dtt (U 2 + sinh 2 ([7) ) + 2 (2i/ - 1) C/(e) - (2z/ - l) 2 log ^~ 



• (B.2) 



Here the dot stands for the t-derivative. Varying (1B.2|) with respect to v yields equation 
U (u) = Since = 0, one has 



An explicit form of Uq{v) is well known [6], 



dx t/o(x) 



f/o = log 2 



a _ 2 , r(i-iz) 



(B.3) 



(B.4) 



Combining (lB~T|) . (IB~3|) with (fB~4|) one arrives to (T2T35|) . 



10 For r _ Q 5 equation (|A.17|) (supplemented by ([OTjl . ([2~28]) . ([235]) ) follows immediately from the results of 
Ref. [7]. It was used in Ref. [13] for the case I = and 2a = 2, 3 . . . . 
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C Appendix: Formula for (S) 



<9£ / r,k 

This section presents some supporting evidences for relation (I3.54p which is expected to hold 
for any £ > 0. 

Using Eqs. fl3.49p . (I3.50p and (I3.53P it is easy to show that for £ > 1 (i.e., in the case without 
soliton-antisoliton bound states) the large-r behaviour of 2) is given by 

2)(r) = -\ cos(2ttA;) K Q {r) + 0(e~ 2r ) , (C.l) 



whereas 



oo smn( 
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CM = ~~\ rVr^ cosh I ^ - 2vrifcj ^(r) + 0(e~ 2 ') . (C.2) 

7r cosh(^) smh(-|-) V 2 / 

Here K u (r) is the modified Bessel function and the symbol 0(e~ 2r ) stands for an asymptotic be- 
haviour of the form oc r 7 e~ 2r as r — > oo with some exponent 7. The leading large-r asymptotic 
of 2) does not depend on £ and hence the L.H.S. of (13. 54f) decays as 0(e~ 2r ). Thus Eq. (13.541) is 
qualitatively consistent with asymptotics 11C.1I) and (1C.2I) . A quantitative comparison can be 
made for £ = 2 and k — ±7. In this case [7] 



(f),,= 2? /"<" t log(i)Aj W+0 (.- 



(C.3) 



(This formula follows immediately from (II. 9p and the large- 1 asymptotic U(t) = - Ko(t) + 
0(e~ 2t )), so that Eq. (l3.54p implies the following easily established identity 



dt t log - f^(t) = / dcu . (C.4) 

W 0W ./.^ 2sinh(7rw) v ; 

Another piece of evidence supporting (I3.54p comes from the consideration of r — > limit. 
Using Eq. (13.351) one finds 

1 /:"-^^^c(-.) + ^c'(o)^(«) + ^, ? ) + o(i), to 

where F(/c,£) stands for 

«-i£7iSi^' (a6) 



2 

.(eft). 



As follows from results of Refs. [26j|27], Q° (a;) is a meromorphic function of the complex 

■|. Also, the product C& (^0 Cfc 



variable k, analytic in the half plane Re(k) > — \. Also, the product C^i 1 ^) (jf^i -00 ) can t> e 



represented in the form 

-(eft)/ w(cft), x_ ™'"v 2 ; / r -if , ,2 



crv) cr'(-u) = —* 2 — — ^— — + 1 + u z x \ , (c.7) 

Sfc v ^ v ; 24sinh(^)cosh(^) V cu 2 + l ) 

32 



where X admits the k — > oo asymptotic expansion as 3i.e(k) > — 



-2n 



(C. 



n=l 



with p = 2k£ and /i = £(1 + £). The coefficient X n (w 2 , /i) are calculated systematically with the 
WKB method applied to the Schrodinger operator (13.371) . They turn out to be polynomials in 
the variables u 2 and h of orders n — 1 and n, respectively. For example, 



X x {u\h) 

X 2 (u 2 ,h) = ^ (J*h'-9h 2 -f h-l) . 
Thus F(k,£) flC.6j) can be written in the form 



(C.9) 



1 + 2£ 



+ 2i log (iT? 



(CIO) 



where ^(z) = ^- logT(z) and F admits the large-fc asymptotic expansion in the half-plane 

&e(Jfe) > -| 



-2;i 



(C.ll) 



n=l 



I have calculated explicitly the expansion coefficients F n up to the order n = 8 and found full 
agreement with the formula 



( d l + 2£ 0_n 



f iog(e(i+o) 



(C.12) 



dt 



sinh(t) e 



-2pt 



1 - 2pt c cff _ 2t 
e 



o t V4tsinh(^)sinh(t(l + 0) 4^(1 + 12 



(ceff = 1 - which follows from Eqs. (13331 . (EE5H , (Q combined with (IZ^IJ) . 



D Appendix: Small-r expansion of c(r, z/ 

Here we explain how to develop the expansion ( I4.38P systematically. 

The partial resummation of the double series (I4.32p yields the following structure 



e2U(t) (8u) 2 X ^ P 3n (X\v) (V 



n=0 

such that 



y P, n {X | U) ( t_yn 

^ X n (l-X) n \8J ' { } 



t 2 (i-xy ^ X n(i-x) 



where X = (j^) 2 t Au , P$ = 1 and P^ n are polynomials in X with the degree deg(P3 n ) < 3n, 



P 3n (X | !/) = (-If Pgn^- 1 | - u) . (D.2) 



33 



With sufficient computer resources, these polynomials can be calculated order by order with 
reasonable facility. Explicitly, 

Ps(X) = y2 2 _ ((l-v)X-v-l)((l-v)X 2 + (Av 2 -2)X + u + l) . (D.3) 



As the next step, one should resummate the series (ID. II) and bring it to the form 



e 2U(t) = t -2 



(D.4) 



where B = B[t) and Vj = Vj(t) , j = —2, —1, . . . are formal power series in t 4 . Note that 
this form is suggested by the Laurent series (I4.33p . Comparing the singular parts at t — r of 
(ILTil) ane (14T331) one finds that 



k = 8u r~ 2v B{r) , (D.5) 

whereas V-2 and V-i are certain differential polynomials of B = B(t): 

V-2 = 16 B 2 (2uB-tB) 2 , (D.6) 
y_ x = 8 (8 v 2 B 2 + (1-8 v)tB B + t 2 B 2 + t 2 BB) . 

Here the dots stand for derivatives with respect to t. To determine the expansion coefficients 
of the formal power series B = B(t) and Vq = Vo(t) one should re-expand (ID.4|) in the powers 
of (X — 1) and compare the terms oc (X — l)~ n for n > with the similar terms from (ID . 1 [) . 
Explicit calculations yield Eq.f l4.38]) and the similar expansion for logK. I verified that the 
power series for c and log k obey the relation (I4.34p up to in twelfth order in r. 
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